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Abstract 

Since the introduction of binomial state as an intermediate state, different intermediate states 
have been proposed. Different nonclassical effects have also been reported in these intermediate 
states. But till now higher order antibunching or higher order subpoissonian photon statistics 
is predicted only in one type of intermediate state, namely shadowed negative binomial state. 
Recently we have shown the existence of higher order antibunching in some simple nonlinear 
optical processes to establish that higher order antibunching is not a rare phenomenon (j. Phys. 
B 39 (2006) 1137). To establish our earlier claim further, here we have shown that the higher 
order antibunching can be seen in different intermediate states, such as binomial state, reciprocal 
binomial state, hypergeometric state, generalized binomial state, negative binomial state and 
photon added coherent state. We have studied the possibility of observing the higher order 
subpoissonian photon statistics in different limits of intermediate states. The effect of different 
control parameters have also been studied in this connection and it has been shown that the 
depth of nonclassicality can be tuned by controlling various physical parameters. 



1 Introduction 

An intermediate state is a quantum state which reduces to two or more distinguishably different states 
(normally, distinguishable in terms of photon number distribution) in different limits. In 1985, such a state 
was first time introduced by Stoler et al. jl]. To be precise, they introduced Binomial state (BS) as a state 
which is intermediate between the most nonclassical number state \n) and the most classical coherent state 
I a). They defined BS as 

M 



\P. M) = Ef=o 5f I") = E \j'''CnP^'{^-pY^i--\n) < p < 1. (1) 

This stat^ is called intermediate state as it reduces to number state in the limit p —> Q and p ^ 1 (as 
|0, M) = |0) and |1, M) = \M)) and in the limit of M ^ oo, p ^ 1, where a is a real constant, it reduces to 
a coherent state with real amplitude. Since the introduction of BS as an intermediate state, it was always 
been of interest to quantum optics, nonlinear optics, atomic physics and molecular physics community. 
Consequently, different properties of binomial states have been studied (2][^. In these studies it has been 
observed that the nonclassical phenomena (such as, antibunching, squeezing and higher order squeezing) 
can be seen in BS. This trend of search for nonclassicality in Binomial state, continued in nineties. In one 
hand, several versions of generalized BS have been proposed |3]|5] and in the other hand, people went beyond 
binomial states and proposed several other form of intermediate states (such as, excited binomial state 
[g], odd excited binomial state j?], hypergeometric state js], negative hypergeometric state |9], reciprocal 
binomial state jio], shadowed state jll], shadowed negative binomial state jl2] and photon added coherent 
state jl3] etc.). The studies in the nineties were mainly limited to theoretical predictions but the recent 
developments in the experimental techniques made it possible to verify some of those theoretical predictions. 
For example, we can note that, as early as in 1991 Agarwal and Tara jis] had introduced photon added 
coherent state as 
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(where m is an integer and la) is coherent state) but the experimental generation of the state has happened 
only in recent past when Zavatta, Viciani and Bellini [l4] succeed to produce it in 2004. It is easy to 
observe that ^ represents an intermediate state, since it reduces to coherent state in the limit m ^ and 
to number state in the limit a ^ 0. This state can be viewed as a coherent state in which additional m 
photon are added. The photon number distribution of all the above mentioned states are different but all 
these states belong to a common family of states called intermediate state. It has also been found that most 
of these intermediate states show antibunching, squeezing, higher order squeezing, subpoissonian photon 
statistics etc. but higher order antibunching has been reported only in shadowed negative binomial state 
[12] . Inspired by these observations, many schemes to generate intermediate states have been proposed in 
recent past (T4][l7 , 



The reason behind the study of nonclassical properties of intermediate states lies in the fact that the most 
of the interesting recent developments in quantum optics have arisen through the nonclassical properties of 
the radiation field only. For example, antibunching and squeezing, which do not have any classical analogue 
[TspO , have extensively been studied in last thirty years. But the majority of these studies are focused 



on lowest order nonclassical effects. Higher order extensions of these nonclassical states have only been 
introduced in recent past (2T|24 . Among these higher order nonclassical effects, higher order squeezing has 



already been studied in detail |21 22 25 26 but the higher order antibunching (HOA) is not yet studied 
rigourously. 

The idea of HOA was introduced by Lee in a pioneering paper j23] in 1990, since then it has been 
predicted in two photon coherent state j23], shadowed negative binomial state [12], trio coherent state |27] 
and in the interaction of intense laser beam with an inversion symmetric third order nonlinear medium |28] . 
Prom the fact that in first 15 years after its introduction, HOA was reported only in some particular cases, 
HOA appeared to be a very rare phenomenon. But recently we have shown that the HOA is not a rare 
phenomenon [29] and it can be seen in simple optical processes like six wave mixing process, four wave 
mixing process and second harmonic generation. To establish that further, here we have shown the existence 
of HOA in different intermediate states, namely, binomial state, reciprocal binomial state, photon added 
coherent state, hypergeometric state, Roy-Roy generalized binomial state and negative binomial state. 

The present work is motivated by the recent experimental observation of intermediate state [l4], theo- 
retical observation of possibility of observing HOA in some simple optical systems [29] and the fact that the 
intermediate states, which frequently show different kind of nonclassicality, form a big family of quantum 
state. But till now HOA has been predicted only in one member (Shadowed negative binomial state) of such 
a big family of quantum states [12]. Motivated by these facts the present work aims to study the possibility 
of HOA in all the popularly known intermediate states. The theoretical predictions of the present study 
can be experimentally verified with the help of various intermediate state generation schemes and homodyne 
experiment, since the criteria for HOA appears in terms of factorial moment, which can be measured by using 
homodyne photon counting experiments [30p3 . In the next section we have briefly described the criteria 
of HOA. In section 3 we have shown that the HOA of any arbitrary order can be seen in BS. In section 4, 
Roy-Roy generalized binomial state [i] is studied and existence of HOA is predicted. Calculational details 
and methodology have been shown only in algebraically simple cases which are described in section 3 and 4. 
Section 5 is divided in several subsections and we have followed the same procedure and have studied the 
possibilities of observing HOA in different intermediate states namely, reciprocal binomial state, negative 
binomial state, hypergeometric states and photon added coherent state. One subsection is dedicated for the 
discussion of one particular intermediate state. Finally section 6 is dedicated to conclusions. 



2 Criteria of HOA 

The criterion of HOA is expressed in terms of higher order factorial moments of number operator. There exist 
several criterion for the same which are essentially equivalent. Initially, using the negativity of P function 
flSl, Lee introduced the criterion for HOA as 



2 



where N is the usual number operator, (A^*^*^) = {N{N ~ 1)...{N ~ i + 1)) is the ith factorial moment of 
number operator, () denotes the quantum average, I and m are integers satisfying the conditions 1 < m < I 
and the subscript x denotes a particular mode. Ba An [27] choose m = 1 and reduced the criterion of ^th 
order antibunching to 



N. 



1 < 



iVi'+i))<(iV«)(iV.). 



(4) 



(5) 



Physically, a state which is antibunched in Ith order has to be antibunched in {I — l)th order. Therefore, we 
can further simplify (|5| as 



i+i 



and obtain the condition for I — th order antibunching as 



< 0. 



(6) 



(7) 



This simplified criterion Q coincides exactly with the physical criterion of HOA introduced by Pathak 
and Garica [28] and the criterion of Erenso, Vyas and Singh [34], recently Vogel has reported a class of 
nonclassicality conditions based on higher order factorial moments [35]. All these criteria essentially lead to 
same kind of nonclassicality which belong to the class of strong nonclassicality according to the classification 
scheme of Arvind et al [36]. Here we can note that d{l) = and d{l) > corresponds to higher order 
coherence and higher order bunching (many photon bunching) respectively. Actually, ^a^'a') = ^A^^'^) is 
a measure of the probability of observing / photons of the same mode at a particular point in space time 
coordinate. Therefore the physical meaning of inequalities (|6} is that the probability of detection of single 
photon pulse is greater than that of a two photon in a bunch and that is greater than the probability of 
detection of three photon in a bunch and so on. This is exactly the characteristic that is required in a 
probabilistic single photon source used in quantum cryptography. In other words all the probabilistic single 
photon sources used in quantum cryptography should satisfy the criteria (|7jl of HOA [3^ . 



3 Binomial State 

Binomial state is originally defined as ([TJ, from which it is straight forward to see that 



a\p, M) 



Similarly, we can write. 



Therefore, 



In - 1 > 



Z^n=0 \ (n-l)!(M-n)!-f^ \^ W J V 

Efio' {M^/^'a-p)''-'-'Y \l > {assurmngn- 1 = I) 
[Mp]^\p,M -I). 

a^\p,M) = [M{M -l)p^]^\p,M -2) 



(8) 



> 



a3|p,M) 



[M{M - 1)(M - 2)p3] 2 |p, M - 3) 



[M{M - 1)....(M - I + l)p']5 |p, M-l) 
\M~l,p). 



(9) 



(M-iy.f 



(M,p|at'^(M-;,p| 



M! 

{M-iy} 



(10) 
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and consequently, 



(M,p|n(')|p,M) = {M,p\a^^a}\p,M) = 



Now substituting (111 in equation ^ we obtain the condition for Ith order antibunching as 



M! 

_{M-l-l)\^ 



i + 1 



< 



(11) 



(12) 



(M - 1)(M - 2)....(M -I) <M^ 



(13) 



which is always satisfied for any M > I and both M and ^ are positive (since every term in left is < M). As 
M is the number of photons present in the field and d{l) is a measure of correlation among {I + 1) photons, 
therefore M > (^+ 1) or M > I. Consequently, a binomial state always shows HOA and the highest possible 
order of antibunching that can be seen in a binomial state is equal to M — 1, where M is the number of 



photon present in the field. From (12 I it is straight forward to see that the number state is always higher 
order antibunched and in the other extreme limit (when p ^ 1, Af — > oo and the BS reduces to coherent 
state) d{l) = 0, which is consistent with the physical expectation. 



4 Generalized Binomial State 

We have already mentioned that there are different form of generalized binomial states (3][5]. For the present 
study we have chosen generalized binomial state introduced by Roy and Roy ji]. They have introduced the 
generalized binomial state (GBS) as 

N 

\N, a,P) = Y^ V^{n,N,a,P)\n) (14) 

n—o 

where, 

^ (a + ^ + 2)^ ^W^l ^ 

with a,P > —1, n = 0, 1, N, and 

(a)o = 1 (a)„ = a(a+ l)....(a + n- 1) . (16) 

This intermediate state reduces to vacuum state, number state, coherent state, binomial state and negative 
binomial state in different limits of a, f3 and N. In order to obtain an analytic expression of d{l) for this 
particular generalized binomial state we need to prove following useful identity: 
Identity 1: 

a(a + 1)„ = (a)„+i (17) 



Proof: Using (16 1 we can write 

/ , IX / , 1^ / , ^ a{a+l)...{a + n) (a)„+i 
(a + Ijn = (a + l)...(a + n) = = — - — , 

Therefore, 



a{a + 1)„ = (a) 



n+l- 



Now it is easy to see that the above identity (171 yields the following useful relations: 

{a+l),+i^{a+l){a + 2)i (18) 

and 

{a + f3 + 2)N = {a + l3 + 2){a + P + 3)n-i = {a + f3 + 2){a + 2 + /3 + 1)n_i. (19) 
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Using (14 1 and (15 I we can obtain 



jV! 

(a+/5+2)jv (n-1)! (N-i 



n-l) 



^N-l r jV(jV-l)! + i (;3+l)jv_i_i \ ^ ,« 



where n = / — 1 has been used. Now we can apply (181 and (19 1 on (20 1 to obtain 



(iV,a,/3|ata|iV,a,/3) = {^^^^ 



(20) 



n\N R\ - / JV(a + l) \ ^ v^AT-l / (Af-1) !(a+2), i \ ^ in 

U|iv,u!,^/ — y(a+l3+2) I ^1=0 \ (q+2+^+1)jv-i/!(A'-1-;)! J ' ' 

where dummy variable I is replaced by n. Therefore, 



(21) 



and 



{N,a,/3\a^'a^\N,a,P) = 



[N{N-1) (jV-;+l)][(a + l)(a+2) ja+l)] 

(a+/3+2)(a+/3+3) (a+/3+i + l) 

_ N\{a+i)\(a+p+l)\ 
~' (Af-/)!Q!(a+/3+; + l)! 



^GBS 



(0 



JV!(a+;+l)!(a+/3+l)! 
(Af-/-l)!a!(a+/3+i+2)! 



f JV(a+l) 1 



i + 1 



J+/3+2) / 

^ [jV(jV~l) (jV-0][(a+l)(a+2) (a+i + 1)] 

(a+/3+2)(a+/3+3) (a+/3+i+2) 



f jV(a+l) \ 
\ (a+/3+2) / 



i+1 



(22) 



The physical condition iV > ? + 1 ensures that all the terms in dilb) are positive. The expression of d(V) is 
quite complex and it depends on various parameters (e.g. a, /3 and TV). Fig 1 shows that for particular 
values of these parameters HOA is possible. As it is expected from the earlier works on the properties of 
HOA [28], the depth of nonclassicality is more in case of dGBs(9) than in dcBsi^)- This is consistent with 
earlier observation. A systematic study reveals that the probability of observing HOA increases with the 
increase of a but it decreases (i.e the probability of higher order bunching increases) with the increase of f3. 
This can be seen clearly in Fig2 and Fig3. Further it is observed (from Figl and Fig2) that for lower values 
of a the probability of bunching increases with the increase of N but for a comparatively large values of a 
(larger compare to /3) the probability of HOA increases with the increase in N [see Fig 2] but the situation 
is just opposite in the case of (3 (as it is seen from Fig3). 
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Figure 1: Higher order antibunching can be seen in Generalized binomial state. Existence of 8th and 9th 
order antibunching (for a = 2 and (3 = l)and variation of depth of nonclassicality with N has been shown. 

While studying different limiting cases of Roy and Roy generalized binomial state, we have observed that 
binomial state and number state always show HOA and d{l) = for coherent state. This is consistent with 
the physical expectation and the conclusions of the last section. 
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Figure 3: Variation of dcssi'^) with f3 and iV for Q!=10. 



5 Other Intermediate States 

As it is mentioned in the earher sections, there exist several different intermediate states. For the systematic 
study of possibiHty of observing HOA in intermediate states, we have studied all the well known intermediate 
states. Since the procedure followed for the study of different states are similar, mathematical detail has not 
been shown in the subsections below. But from the expression of d{l) and the corresponding plots it would 
be easy to see that the HOA can be observed in all the intermediate states studied below. 

5.1 Reciprocal binomial state 

Reciprocal binomial state (RBS) can be defined as jio] 

= 1 E (k) <^''^''^^'Hk) (23) 

fc=0 ^ ^ 



where A'^ is a normalization constant. Using antinormal ordering and procedure adapted in section 3, we can 
obtain; 

^ fi\ - \^ ( 1^i (i+i)!^ {N+{i+i-iy. 



1=0 



(24) 



TT csc(TT{l+N))(l + l)\^Gamma(l + l-N} _ at' + I 
Mt{Gamma{2+l}Gamma{-Nyf 



Where Gamma denotes the Gamma function. The possibility of observing HOA in reciprocal binomial 
state can be clearly seen from the Figure 4. But it is interesting to note that the nature of singularity and 
zeroes present in the simplified expression of duBs{l) as expressed in the last line of (24 1 can provide us 
some important information. For example, the underlying mathematical structure of the criterion of HOA 
and that of reciprocal binomial state demands that I and N be integers but if both of them are integer 
then djiBsiO has a singularity as the csc[7r(Z + N)] term present in the numerator blows up. But this local 
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singularity can be circumvented by assuming I integer and N integer. In this situation (i.e. when I 
and N tends to integer value) (A^*^'^) is finite and consequently d{l) is also finite. This is the reason that the 
singular nature of the simplified expression of djiBs{l) is not refiected in the Figure 4. In the analysis of the 
d-RBsil) it is also interesting to observe that Gamma{—N) = oo for N — integer^ and in an approximated 
situation when I — > integer and M = integer, the csc[7r(Z + N)] term in the numerator is no more singular 
and as a result(A^(')) = and diiBs = —N'^^^^y In this situation one can observe HOA for arbitrarily large 
values of I and N. Thus physically, it is expected that in reciprocal binomial state higher order antibunching 
of any arbitrary order will be seen and HOA will not be destroyed with the increase of iV, as it happens (for 
some particular values of a and (3) in the case of generalized binomial state (see Figure 1). 
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Figure 4: Variation of duBsi^) and dj^Bsi^) with photon number N 



5.2 Negative Binomial State 

Following Barnett |38] we can define Negative Binomial state (NBS) as 



where Cn{ri,M) 



\7j,M)=J2Cn{v,M)\n) (25) 

n=M 
-,1/2 



< ?7 < 1 and M is a non-negative integer. This 



n 
M 

intermediate state interpolates between number state and geometric state. Following the mathematical 
techniques adopted in the earlier sections we obtain 



dNBs{l) = V j^, j (26) 

where 2Fi{a,h\c\ z) is a conventional hypergeometric function. Variation of d{l) with various parameters 
such as 77, I and M have been studied and are shown in Figure 5- Figure6. From these figures one can observe 
that the state is not always antibunched, rather the plot of d^Bsi^) has a very sharp rise near 77 « .15 and 
M « 10. From Figure 5 we can observe that the possibility of higher order bunching get destroyed with the 
increase of M or rj. The second feature has clearly been shown in Figure 6, where we can see that dNBsi^) 
become positive for a small region for M — 10 but for a bigger region (77 « .3 to 77 = 1) it remains negative 
and thus shows the existence of HOA. The broad features remain same for the other orders (other values of 
I) of antibunching. 

In the limit iV ^ the negative binomial state reduces to geometric state (GS). In this limit d^Bsil) 
reduces to 

dG5(0-^((l-'/)'+Sa+l)!-l)- (27) 

It is interesting to observe that the above expression has a singularity at 7; = and dasil) —00 in the limit 
7/ 0. Consequently one negative values of das can be seen at very very small values rj but this is not the 
signatures of HOA, rather this is the signature of the existence of a strong singularity in the neighborhood. 
This can further justified by the fact that for any finite value of I there does not exist any real root, (whose 
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Figure 5: Variation of dNBsi^) with 77 and M 





Figure 6: Variation of dNBsi^) with rj for M — 10 

value is close to zero or which is negligibly small compared to 1) of dcsil) — 0. Thus there is no oscillation 
between bunching and antibunching. We further observe that c?gs(0 ~^ ~1 in the limit 77 ^ 1 and for I > 3 
there exists only one physically acceptable real root of das{l) = 0. By physically acceptable real root we 
mean that it lies in [0, 1]. Before this value of 77 (or before the physically acceptable real root) the state 
shows higher order super poissonian photon statistics but immediately after the root it becomes negative 
and thus shows HOA or higher order subpoissonian photon statistics. As we increase I the real root shifts 
in right side of the real axis (i.e towards 77 = 1). It can be clearly seen in the Figure 7. From this figure it 
can be easily seen that it satisfies all the physical properties of HOA derived in (281. 



5.3 Photon added coherent state 

Photon added coherent state ([SJ or PACS, which was introduced by Agarwal and Tara [l3] can be defined 

as 

\a\^ \ 00 n 



I \ _ exp(-LJ-) 

v/£™(-|aP)m!^„ n\ 
where is Lauguere polynomial of mth order. Rigorous operator algebra yields 

cxp ( — a 



(28) 



dpAcs{^) 



(m!(( + l)!)2 iFi(-m;l 

cxp(-Q^)(-m+rn iFi (l+m; ) + (l+m)a^ iFi(2+m;2;a^)) 
lFi(-m;l;-a2) 



(29) 



where, pFq is the generalized Hypergeometric function. The analytic expression for dpAcs(0 quite 
complicated and it is difficult to conclude anything regarding its photon statistics directly from ( [29| but 
we have investigated the variation of dpAcs with a, I and m and could not find any region which does not 
show HOA. Therefore, HOA can be seen in this particular intermediate state. This fact is manifested in Fig 
8 and Fig 9. From these figures it is easy to observe that depth of nonclassicality increases monotonically 
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Figure 7: Variation of (iGs(8), dQs{Q) and rfGs(lO) with respect to -q. 

with the increase of m and L The variation of depth of nonclassicality with a has a deep for a small 
value of a (see Fig 8. and Fig. 9). Although dpAcs is always negative, initially its magnitude {dpAcs{l) 
without the negative sign) increases, then decreases and then become a monotonically increasing function. 
Actually for the smaller values of a, an effective contribution from the combination of all the hypergeometric 
functions appears and dominates but as soon as a increases a bit, the exp(— a^) term starts dominating 
and as a consequence depth of nonclassicality increases monotonically. Here we would also like to note 
that in contrast to the photon added coherent state, d{l) is always positive for the analogous state |a,— m), 
introduced by Sivakumar |39]. Thus |a, — m) always shows higher order superpoisonian photon statistics. 
Further, we would like to note that photon added coherent state which, is intermediate between coherent 
and fock state has already been experimentally generated in 2004 jl4]. Therefore, it is technically feasible 
to observe higher order antibunching for an intermediate state. 




Figure 8: Variation of dpAcsi^^) with a and m 



5.4 Hypergeometric state 

Following js] hypergeometric state (HS) can be defined as 

M 

\L,M,rj) = Y.^n{V:L)\n) (30) 

n=0 
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Figure 9: Variation of 10dpACs(3) and dpACs(4) with a for m = 15, the soHd line denotes (ipACs(3) and 
the dashed Hne denotes rfp^c's(4) to keep show the variation in the same scale dpACS^^ is multiplied by 
10. The plot shows that depth of nonclassicality of dpACs(4) is always greater than that of dp^cs(3) which 
is consistent with the properties of HOA 



where the probability jj is a real parameter satisfying 0<?7<l.Lisa real number satisfying 

L > max {Mr]~^, M(l - jy)"^} , 



Lr] 
n 



M -n 



L 

M 



and 



a{a — l)...(a — n + 1) 



n! 



1. 



(31) 
(32) 

(33) 



Here a is not necessarily an integer. Using the techniques adopted in the earlier sections and a bit of operator 
algebra we can obtain a closed form analytic expression for d{l) as 



dHs{l) = - {Mriy+^ + 



{L - I - iy.M\{Lri)l 



L\{M- 



l)l{Lr) - 1 - l)\ 



(34) 



From Fig. 10 it is clear that HOA can be observed in hypergeometric state. It is also observed that the 
depth of nonclassicality increases with the increase in rj and M. Hypergeometric state reduces to binomial 
state, coherent state, number state and vacuum state in different limits of M, L and rj. It has been verified 
that if we impose those limits on dns then we obtain corresponding photon statistics. 



6 Conclusions 

In essence all the intermediate states studied in the present paper show HOA i.e. higher order subpoissonian 
photon statistics. But it neither mean that all the intermediate states are higher order antibunched (for 
example \a, —m.) is always higher order bimriied) nor an intermediate state which shows HOA will show it 
for all possible values of the control parameters (for example negative binomial state and generalized binomial 
state shows both higher order bunching and higher order antibunching for various parametric values). Thus 
we can conclude that, as far as HOA is concerned there does not exist any common characteristics among the 
different intermediate states but most of them show HOA. Further, we have seen from (Figurel-FigurelO) 
that the depth of nonclassicality of a higher order antibunched state varies with different control parameters 
(e.g. a, N, m etc.). These parameters represent some physical quantity and their value may be controlled and 
consequently by controlling these parameters we can control the depth of nonclassicality. Photon statistics 
(factorial moment) of an intermediate state can be obtained experimentally by using homodyne detection 
(photon counting) technique. These facts along with the recent success in experimental production of 
intermediate state open up the possibility of experimental observation of HOA in intermediate state. Thus 
the present work strongly establishes the fact that HOA is not a rare phenomenon. 

The prescription followed in the present work is easy and straight forward and it can be used to study 
the possibilities of observing higher order antibunching in other intermediate states (such as negative hyper 
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Figure 10: Variation of dnsi^) with rj and M, when lowest allowed values of L have been chosen at every 
point. 

geometric state, excited binomial state and odd excited binomial state) and other physical systems. Thus 
it opens up the possibility of studying higher order nonclassical effects from a new perspective. This is 
also important from the application point of view because any probabilistic single photon source used for 
quantum cryptography has to satisfy the condition for higher order antibunching. 

Acknowledgement: AP thanks to DST, India for partial financial support through the project grant 
SR\FTP\PS-13\2004. 
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